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ABSTRACT 



Context. Molecular outflows from high-mass young stellar objects provide an excellent way to study the star formation process, and 
investigate if they are scaled-up versions of their low-mass counterparts. 

Aims. We selected the nearby massive star forming region IRAS 17233-3606 in order to study the kinematics and physics along the 
molecular outflow(s) originating from this source. 

Methods. We observed IRAS 17233-3606 in CO, a typical tracer of gas associated with molecular outflow, with the Submillimeter 
Array in the (2-1) transition, and with the APEX telescope in the higher excitation (6-5) line. Additional infrared H2 observations 
were performed with the UKIRT telescope. The CO data were analysed using a LVG approach. 

Results. Our data resolve the previously detected molecular outflow in at least three different components, one of them with a high 
collimation factor (~ 4), and characterised by emission at extremely high velocities (|v - v LSR | > 120 km s -1 ), The estimate of the 
kinematical outflow parameters are typical of massive YSOs, and in agreement with the measured bolometric luminosity of the source. 
The kinematic ages of the flows are in the range 10 2 - 10 3 yr, and therefore point to young objects that still did not reach the main 
sequence. 

Key words. ISM: jets and outflows - ISM: molecules - Stars: individual: IRAS 17233-3606 - Stars: formation 



1. Introduction 

Although massive stars are fundamental for astrophysics, little 
is known of their earliest evolutionary phases. The main prob- 
lem for the understanding of massive star formation comes from 
the radiation pressure that massive stars exert on their surround- 
ing medium once they have reached a mass of around 8 Mq, 
that should prevent further accretion onto the star. As a conse- 
quence, stars with mass larger than 8 Mq should not exist ( |Palla| 



ity among early spectral type stars (e.g., [Preibisch et aL 2001) 



& Stahler T993J, in clear disagreement with observations. Two 



make observations of these phases of massive YSOs a challenge, 
especially at (sub)mm wavelengths where current facilities still 
lack the necessary resolution to resolve single massive objects at 
such distances. 

Despite the observational effort in the last decade to in- 
vestigate the early evolutionary phases of massive stars, only 
few sources show typical phenomena associated with low-mass 
star formation such as collimated jets and outflows, circumstel- 



main theoretical scenarios have been proposed to solve this prob- 
lem: (1) accretion as in the low-mass case through more mas- 
sive circumstellar disks and/or sufficiently high accretion rates 
(e.g., |Yorke & Sonnhalter| |2002| l and (2) coalescence of lower- 
mass stars belonging to the cluster that harbours the high-mass 
star(s) HBonnell et aL| [19981 |Bonnell & Bate] |2002> . One of the 
fundamental observational tests to distinguish between the dif- 
ferent proposed models is to detect collimated outflows and ac- 
cretion disk from massive young stellar objects (YSOs), and to 
assess the fraction of still embedded massive YSOs with disks 
( |Cesaroni et"aT] |2007) |Krumholz & Bonnell||2007] l. Since mas- 
sive pre-main sequence stars may loose any disk they originally 
had due to the radiation field (e.g., |Natta et aT| |2000} [Fuente 



lar disks even at high angular resolution (Beuther & Shepherd 
|2Tj05l|Arce et al.||2TjO"7l|Cesaroni etaL||2007| l. In addition, nei 
ther collimated jet-like outflows or circumstellar disks have been 
detected in sources more massive than early B- and late O-type 
YSOs, although recently Zapata et al. ( 2009 1 detected a rel- 



atively collimated outflow (14°) towards W51 IRS2 and sug- 
gested the presence of a Keplerian infalling ring around a central 
object of at least 60 M Q . 

It is still not clear whether the lack of detection of collimated 
outflows and circumstellar disks around O-type YSOs is an ob- 
servational bias due to the poor angular resolution and sensitivity 
of the observations or whether it reflects a real difference in the 



way massive stars form. Beuther & Shepherd (2005) suggested 



et al. 2003 1, such studies are best done in earlier evolutionary 



phases. However, large distances and a high degree of multiplic- 
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an evolutionary scenario in which the collimation of the out- 
flows decreases with time, due to the interaction with the wind 
from the central object. Hence, the need to investigate high-mass 
YSOs in different evolutionary stages and of different spectral 
types at high linear resolution is clear. 
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The region harbouring the prominent far-infrared source 
IRAS 17233-3606 (hereafter IRAS 17233) is one of the best 
laboratories in which to attack the problem of massive star for- 
mation. IRAS 17233 first came to attention through its very in- 
tense H 2 Q, OH, and CH 3 OH m asers (e.g., |Caswell et aL][T^80t 

TWT) . The lumi nosity of IRAS 17233 2. Observations 



(e.g., |Beuther et ah] |2002a| |Leurini et al.| |2006[ |van Kempen| 
et al.| |2009[ l because of the low transmission of the atmosphere 
at these frequencies. 



Fix et al. 1982 Menten 



range s between 1.4 x 10 3 L Q (d = 2.2 kpc, Hughes & MacLeod 
1993| based on the IR AS colours) and 1.7x 10 4 L Q (d = 0.8 kpc, 
Faundez et al.| 2004| obtained by integrating the spectral energy 
distribution). Previous studies seem to agree that IRAS 17233 
is located at the near kinematic distance (between 700 pc and 



2.2 kpc, |Miettinen et al. 



2006] |Forster & Caswell 
~ 16 kpc, e.g 



T989> rathe r 
MacLeod et al.| | 1998) . 



than at the far distance ( 

This is suggested by the high measured intensities of continuum 
and line emission at practically all the observed wavelengths, 
which would indicate exceedingly high luminosities if the source 
were at the far distance and, persuasively, by the fact that it is at 
an angular distance of more than 0.5 degrees below the Galactic 
plane. Moreover, the source lies at the edge of an infrared dark 
cloud (Leurini et al., in prep.). In the following analysis, we as- 
sumed a distance of 1 kpc for IRAS 17233. 

Recently, high angular resolution multi-radio wavelengths 
observations by Zapata et al. (2008) resolved a cluster of nine 
compact radio sources of different nature in the region. Of these 
objects, one corresponds to the Hn region already mapped at 
cm-wavelengths by different authors (e.g., Hughes & MacLeod] 
|1993[|Walshet al. 1998 ), while four others are found close to the 
maser zone (i.e., the area of the region where the H2O, CH3OH 
and OH masers are detected, see Fig. [TJ. Zapa ta et al.| ((2008 ) 
also found that one of these four sources (labelled VLA 2d in 
their paper) is at the centre of a bipolar north-south outflow de- 
tected by |Fish et aL] ( |2005) in OH masers. [Zapata et al.|p008) 
also reported the presence of a ring-like structure traced by H2O 
masers, which they interpret as due to northeast-southwest out- 
flows. Leurini et al.| ( |2008) reported the discovery of a bipolar 
outflow originating from the vicinity of IRAS 17233, through 
low resolution observations performed in CO(3-2) with the 
APEX telescope (~ 18" resolution). From these data however, 
the powering source of the molecular outflow remains unidenti- 
fied. [LeunnTetal. (2008 ) also reported a rich molecular spectrum 
from a position very close to the centre of the outflow, typical 
of hot molecular cores near massive young stellar objects, that 
represent the stage of massive star formation when the newly 
formed star is efficiently heating the surrounding medium up to 
hundreds of K, but it has not yet developed an ionised region 
(e.g. 



Kurtz etal. 2000 Cesaroni et al. 2003 1. 



In this paper, we present observations of the molecular out- 
flow in IRAS 17233 at 220 and 230 GHz in 13 CO(2-l) and 
CO(2-l) with the Submillimeter Array (SMA), and at 690 GHz 
in CO(6-5) with the APEX telescopaH The corresponding con- 
tinuum data of the SMA observations, as well as the other spec- 
tral line data, will be presented in a forthcoming paper. While 
the SMA CO(2-l) data allow us to reach a high angular resolu- 
tion (see Table [T| and study the outflow in its detail, the APEX 
CO(6-5) observations trace the warm component of the outflow, 
and allow a reliable estimate of its excitation conditions. Despite 
the importance of observations of high rotational CO transitions 
in the study of molecular outflows, only a handful of sources 
have been mapped in rotational CO lines with J higher than 3 



1 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



2.1. SMA observations 

We observed IRAS 17233 with the SMA interferometer on April 
10, 2007 in the compact configuration with seven antennas. The 
receivers operated in a double-sideband mode with an IF band 
of 4-6 GHz so that the upper and lower sideband were sepa- 
rated by 10 GHz. The central frequencies of the upper and lower 
sideband were 220.4 and 230.4 GHz respectively, allowing to 
detect the I3 CO(2-l), C lg O(2-l) and SO(5 6 - 4 5 ) lines in the 
lower side band of the receivers, the CO(2-l) transition in the 
upper side band. The correlator had a bandwidth of ~ 1.9 GHz 
and the channel spacing was 0.406 MHz, corresponding at this 
frequency to a velocity resolution of 0.5 km s _1 . The adopted 
systemic velocity was -3.4 km s . The observations were per- 
formed under good weather conditions with zenith opacities 
t(225 GHz) between 0.1 and 0.14 measured by the National 
Radio Astronomy Observatory (NRAO) tipping radiometer op- 
erated by the Caltech Submillimeter Observatory (CSO). 

We covered the molecular outflow from IRAS 17233 pre- 
sented by Leurini et al. ( |2008) with a mosaic of two 
fields centred at a 2 ooo=17"26 m 41*.76, <5 20 oo=-36°09'0'.'5 and 
a 2 ooo=l7 h 26 m 42 s .4, (5 2 ooo=-36 09'ia.'5. The projected base- 
lines ranged between 5 and 88 m. The short baseline cutoff im- 
plies that source structures > 54" are filtered out by the obser- 
vations. However, intermediate structures are also affected by 
missing fluxes. 

Bandpass calibration was done with 3C454.3. We used 
Callisto for the flux calibration which is estimated to be accurate 
within 20%. Gain calibration was done via frequent observations 
of the quasars 1626-298 and 1713-269. Measured double-side 
band system temperatures corrected to the top of the atmosphere 
were between 100 and 400 K, mainly depending on the eleva- 
tion of the source. Details on the observational setup are given 
in Tabled] 

The initial flagging and calibration was done with the IDL 
superset MIR originally developed for the Owens Valley Radio 
Observatory and adapted for the SM Ap] The imagi ng and data 

" 7995) and 
weighting. 
1] The high 



analysis was conducted in MIRIAD (Sault et al. 
MAPPING^ Images were produced using natura 
The resulting synthesised beams are listed in Table 
degree of elongation of the beam is due to the low elevation of 
the source from Hawaii. However, the spatial resolution is high- 
est along the direction perpendicular to the outflow, making the 
dataset well-suited to study its morphology and kinematics. In 
order to compare the SMA observations to the APEX data, we 
resampled the final data to a spectral resolution of 2 km s ; the 
corresponding r.m.s. is 0.05 Jy/beam. 

2.2. APEX observations 

IRAS 17233 was observed in the CO(6-5) transition with the 
CHAMP + dual-frequency heterodyne submillimetre array re- 
ceiver on APEX. The observations were performed on 2007, 
October 20 and 22, under good weather conditions (with 0.9 



2 The MIR cookbook by Charlie Qi can be found at http://cfa- 
www.harvard.edu/~cqi/mircook.html 



J http://www.iram.fr/IRAMFR/GILDAS 
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Fig. 1. Overview of the region: the black contours show the 1.3 cm continuum emission of IRAS 17233 (Zapata et al. 2008 1, from 
0.12 Jy beam 1 (3cr) in step of 0.24 Jy beam -1 . The Hn region and the maser zone are labelled. Overlaid on the map are the 
positions of the OH (blue and red squares, [Fish et al. 2005 I, CH3OH (black triangles, [Walsh et al. 1998) and H2O (blue and red 



open triangles, Zapata et al. 



as reference for the APEX observations 



2008) masers. The cross marks the position of the peak of the SMA 1.3 mm continuum emission, used 



Table 1. Observational parameters. 



Transition Frequency HPBW P.A. Av 

(GHz) O (°) (kms- 1 ) 

lj CO(2-l) 220.3987 5.43x1.88 30 050~ 

CO(2-l) 230.5380 5.45x2.08 28 0.50 
CO(6-5) 691.4731 8.9 - 0.05 



is not uniform over the channels, because the northern part of the 
map was observed only on October 20, and because of the dif- 
ferent weather conditions during the two days. The r.m.s. noise 
is of the order \-A K, larger on the northern part of the map. 
The maps were produced with the XY_MAP task of CLASS90, 
which convolves the data with a Gaussian of one third of the 
beam: the final angular resolution of the data is 974. 



and 0.7 mm precipitable water vapour). The peak of the 
SMA continuum emission (ffnooo = 17 /, 26'"42.455 J , feooo = 
-36°09' 18.047") was used as reference position. The system 
temperature was between 4500 and 14000 K on October 20th, 
between 2800 and 4500 K on October 22nd. Pointing was estab- 
lished by total power continuum cross-scans on SgrB2(N), and 
found to be accurate within ~ 3". 

The beam efficiencies of CHAMP + were determined via ob- 
servations of Mars and Jupiter, the first at 691 GHz, the second 
at 661 GHz, and found to be 0.38 and 0.45, respectively. We 
therefore adopted a beam efficiency of 0.40 to convert from an- 
tenna temperature into main-beam temperature units. The origi- 
nal spectral resolution of the data (0.05 km s _1 ) was smoothed to 
2 km s _1 for our analysis. The r.m.s. noise in the final data cube 



2.3. UKIRT observations 

We obtained near-infrared wide field images using WFCAM on 
UKIRT through JHK broad band and a H2 narrow band filter on 
May 3rd and 29, 2008. Integration times were 3.75 min for the 
broad band filters and 24 min for the narrow band filter. We used 
the data as reduced through the CASU pipeline. We used the K 
broad band image to discriminate between continuum and emis- 
sion line features in the H2 narrow band image. We created a 
continuum subtracted, pure H2 emission line image by properly 
registering the broad and narrow band images and smoothing 
the narrow band image to match the poorer seeing of the broad 
band image. We derived a flux scaling factor by comparing the 
counts of a number of stars in both filters and subtracted the 
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scaled broad band image from the narrow band image. The re- 
sulting emission line image reveals a number of jet-like features, 
of which we will only discuss the counterparts of the CO out- 
flows here. The other flows from the cluster forming core and 
along the filamentary dark cloud will be presented elsewhere 
(Stanke et al., in prep.). 



3. Observational results 

3.1. CO emission 

Our new observations significantly improve the angular resolu- 
tion of the previous study (Leur ini et al.| |2008| l, revealing the 
structure of the molecular outflow in detail. The CO transitions 
are characterised by emission up to extremely high velocities in 
the red- and blue-shifted wings. We a dopt here a systemic ve - 
locity of -3.4 km s _1 for the source ( [Bronfman et al. 1996 1, 
which is confirmed by spectral features in the SMA dataset with 
peak around -3.0 km s . The blue-shifted emission extends to 
~ -200 km s" 1 in the CO(2-l) line, and to ~ -80 km s" 1 in the 
CO(6-5) transition. The red-shifted lobe extends to velocities up 
to 120 km s" 1 in CO(2-l), to ~ 60 km s" 1 in the (6-5) line. In 
Fig. [2] we show the CO spectra towards the peak of the 1.3 mm 
continuum emission. 

Based on the channel maps, we identify two different ve- 
locity regimes in the CO(2-l) emission: the extremely high ve- 
locity regime (EHV: blue: -200 < v L sr < -130 km s _1 ; red: 
90 < vlsr < 120 km s _1 ), and the high velocity regime (HV: 
blue: -130 < v L sr < -25 km s _1 ; red: 16 < v LS r < 50 km s _1 ). 
The integrated intensity in the two velocity intervals is presented 
in the left panel of Fig. [3] For the high velocity red-shifted emis- 
sion, we did not integrate over the velocities (50,90) km s _I 
because of contamination from other molecular species at the 
peak of the mm continuum emission (see Fig. [2]). We stress here 
that the molecular spectrum towards the peak of the continuum 
emission is extremely rich and typical of a hot-core, with very 
few channels free of line emission. For this reason, we cannot 
exclude contamination from other molecular species, but this 
problem should affect our maps only at the centre of the mm 
continuum emission. The negative platforming seen in the spec- 
tra towards this position (Fig. |2]i is probably a result of at least 
two effects: first, it is not unlikely that we overestimate the con- 
tribution of the continuum, since at the position of the hot core 
the density of lines is extremely high and it is difficult to define 
channels free of line emission. Since the negative platforming 
affects only the area of the map where the continuum emission 
is detected (see Fig. |4] and [5] for comparison of spectra at other 
positions), we believe this to be the most likely reason for it. 
However, higher energy transition lines (e.g., CO(6-5) but also 
CO(3-2), Leurini et al.| [2008) have wing emission at very high 
velocities, and therefore it is possible that large structure emis- 
sion is filtered out by the SMA and results in negative features. 

The SMA CO maps reveal a clumpy structure with well sep- 
arated blue- and red-shifted emission, and an overall structure 
well aligned along the N-S direction and centred close to the 
peak of the mm continuum emission (aj2ooo = 17 /, 26™42.455 I , 
5j20oo = -36°09'18'/047, see Fig. [6]>. This peak is associated 
with the maser zone and not with the Hn region reported by sev- 
eral authors (see Fig. [TJ. The red-shifted emission shows four 
peaks: Rl associated with HV and EHV emission, R2 and R3 
only with HV gas, and R4 with EHV emission only. Similarly, 
the blue-shifted emission peaks at four positions: Bl and B2, as- 
sociated only with HV gas, and B3 and B4 associated with EHV 
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Fig. 2. Spectra of CO and its isotopologues at the peak of the 
mm continuum emission; the spectral features due to different 
molecular species are marked. As discussed in Sect. 3.1, the neg- 
ative platforming, which occurs only in the region where contin- 
uum emission is detected, is probably due to an overestimate 
of the continuum. In addition, absorption features are probably 
due to filtering of large structures. Note also that several hot- 
core molecular tracers have been detected. Such emissions will 
be analysed and discussed in a forthcoming paper focused on the 
close surroundings of the YSOs. 



emission. The spectra at the eight peaks of the CO(2-l) emission 
are shown in Figs.|4]|5] 

Figure [3] also shows the integrated emission of the blue- and 
red-shifted wings of the CO(6-5) line. Despite the different res- 
olutions, the two maps are very similar. Emission is detected at 
all positions detected in CO(2-l), except on the high-velocity 
B4. 

To roughly quantify the degree of collimation of the whole 
outflow structure as traced by CO(2-l), we divided the length 
of the outflow by its width, which results in a collimation factor 
f c ~ 4. In doing so, we derived the length (B4-R4) and width 
(B4-B3) of the flow from the 5cr level contour of the EHV in- 
tegrated intensity, and assumed that this comes from only one 
molecular outflow (see § |4] for a discussion on the multiplic- 
ity of flows in the region). Similarly, the opening angle of the 
outflow can be estimated from the length of the deconvolved 
major axis a and the length of the deconvolved minor axis b, 
6 = 2xarctan (bid). The derived value is 6 ~ 11°. Alternatively, 
if the H2 emission in the south of R4 is associated with the same 
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Fig. 3. Integrated emission of the blue- and red-shifted wings in the CO(2-l), and (6-5) lines, and continuum emission at 4.5 //m 
overlaid on the H2 emission (grey scale). For the CO(2-l) transition, the solid contours show the EHV blue- (v=[-200,-130] km s -1 ) 
and red-shifted emission (v=[90,120] km s _1 ); the dashed contours mark the HV blue- (v=[-130,-25] km s _1 ) and red-shifted 
emission (v=[16,50] km s _1 ). For the CO(6-5) line, the solid contours show the blue- (v=[-80,-15] km s ) and red-shifted emission 
(v=[8,60] km s _1 ). All contours start from a 5cr level and are in step of 5cr (see Table|2]for the r.m.s noises). The squares mark the 
positions of the red-shifted, blue-shifted CO(2-l) peaks, the triangle the CH3OH masers. The beam of each observational dataset 
is shown in the bottom left corner. In the left panel, the dashed lines and the solid arrow outline the possible molecular outflows 
discussed in §[4]and labelled IRAS 17233-OF1, -OF2, -OF3. The dotted line marks the possible connection between R3 and B3. 
The solid lines in the left panel outline the region mapped in Fig. [6] 



outflow, the collimation factor would be ~ 6 and the opening an- 
gle 8°. Assuming a distance of 1 kpc, the projected extension of 
the outflow on sky is 0.2-0.3 pc. 



3.2. CO outflow vs. H 2 emission 

Figure [3] compares the red- and blue-shifted CO emission with 
the H2 at 2.12 fim map and the 4.5 /im continuum emission 
from the Spitzer Space Observatory (Benjamin et al. 2003 ). This 
emission is often attributed to H2 lines tracing shocked gas aris- 



ing from outflow activity ( |Noriega-Crespo et al. |2004[ |Ybarra| 
& Lada 2009), although contamination from rovibrational lines 
of CO may also happen in hot gas (Marston et al. 2004| l. H2 
emission at 2.12 fim is found in correspondence to all the blue- 
shifted clumps. Emission at 4.5 /mi is detected towards both B3 
and B4. Interestingly, the EHV B3 and B4 clumps are located at 
the edges of the IR bow-structure traced by H2 and by the 4.5 /mi 
emission, and could trace material recently shocked. 

Probably due to extinction, only Rl among the red-shifted 
clumps shows an H2 counterpart (see discussion in paragraph 
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Fig. 4. Spectra of CO and its isotopologues at the peaks of the CO(2-l) blue-shifted emission. The dashed lines mark the velocities 
-200, -130 and -25 km s used to define the EHV and HV regimes. Absorption features are probably due to filtering of large 
structures. 



Table 2. Noise level in the integrated intensity maps presented 
in Fig. [3] and [7] 



Transition rms 

[Jy beam -1 ] 



CO(2-l) 


red-shift., HV 


0.5 


CO(2-l) 


red-shift., EHV 


1.3 


CO(2-l) 


blue-shift., HV 


8.0 


CO(2-l) 


blue-shift., EHV 


1.3 


CO(6-5) 


red-shift. 


1035 


CO(6-5) 


blue-shift. 


863 


B CO(2-l) 


red-shift. 


0.3 


B CO(2-l) 


blue-shift. 


1.0 


C ls O(2-l) 


red-shift. 


0.4 


C ls O(2-l) 


blue-shift. 


0.7 


SO(5 6 - 4 5 ) 


red-shift. 


0.4 


SO(5 6 - 4 5 ) 


blue-shift. 


2.0 



§|4j. Figure[3]shows also H2 emission ~ 40" south of the YSOs, 
but it is not clear whether it is tracing emission related to the 
red-shifted IRAS 17233 outflow. The EHV Rl clump lies ahead 



of an independent, short H2 jet (see Fig. [6}: also in this case it 
seems we are observing material accelerated along the jet traced 
byH 2 (see§|g. 

3.3. Emission from other molecular species 

Together with the main isotopologue of CO, also the 13 CO, 
C ls O(2-l), and SO(56 - 4s) transitions show non-Gaussian 
wings in the spectrum taken at the position of the mm contin- 
uum emission (Fig.|2|i. Figure|7]shows the integrated intensity of 
the three transitions in the red- and blue-shifted wings. For the 
red-shifted emission 13 CO emission, we only used the range of 
velocities (15,20) km s because of contamination from other 
molecular species. 

The emission traced by the three molecular species is more 
compact than that traced by 12 CO, as expected. For the red- 
shifted emission, the 13 CO emission extends towards the fea- 
tures R2 and R3, while the SO(56 - 4s) transitions has a 5cr de- 
tection on R2 and extends towards the Rl position, where EHV 
red-shifted CO(2-l) and H2 are detected. For the blue-shifted 
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Fig. 5. Spectra of CO and its isotopologues at the peaks of the CO(2-l) red-shifted emission. The dashed lines mark the velocities 
120, 90, 50 and 16 km s _1 used to define the EHV and HV regimes. Absorption features are probably due to filtering of large 
structures. 



emission, all transitions are detected on Bl, although none of 
them peaks at this position. All three transitions show two spa- 
tially separated blue-shifted emission peaks; the same behaviour 
is seen in the CO(2-l) line at velocities between -40 and -20 
km s _1 , similar to the range of velocities used to obtain the blue- 
shifted integrated maps of 13 CO and SO. However, the position 
of the two peaks is not the same for the four species. The east 
blue-shifted peak is probably associated with an outflow traced 
by a H2 jet and by the H2O masers (see discussion in the follow- 
ing paragraph); the west peak is probably related to one of the 
other two outflows of the region. 



4. Multiple molecular outflows from IRAS 17233 

A zoom in the central region of the molecular outflow reveals 
a complex picture. Figure [6] shows the continuum emission at 
1 .3 mm towards the inner region of the outflow detected in CO; 
overlaid on the map are the 1.3 cm continuum emission from 
|Zapata et al. (2008 ), the positions of the OH H 2 Q and CH 3 OH 



Four compact sources are detected in cm continuum emis- 
sion within ~ 2" from the peak of the mm continuum emission 
(and labelled VLA 2a, 2b, 2c and 2d by |Zapata et alT] [2008) . 
The brightest source, VLA 2a, is at C/7 from the mm continuum 
peak, suggesting that the two sources probably coincide. The cm 
VLA 2a emission reveals two components, one dominated by a 
classical Hn region, the other likely due to dust emission. On 
the other hand, VLA 2b shows a negative spectral index, possi- 
bly associated with synchrotron emission. Finally, VLA 2c and 
2d have a spectral energy distribution compatible with a hyper- 
compact Hn region (HCHn) or with a massive dusty core or disk. 

Other information on the sources comes from the maser 
emission. The CH3OH masers (black triangles in Fig.|6|, so far 
detected only towards massive stars (e.g., |Xu et al.||20 08), seem 
to indicate their association with VLA 2b, although with an un- 
certainty of ~ 1" in the absolute position. A bipolar outflow is 
detected in OH, with blue-shifted emission towards the north and 
red-shifted emission towards the south as for the CO emission. 
OH masers (filled squares) are observed from molecular mate- 
rial surrounding ultra-compact Hn regions, and from protostellar 



maser spots ( Fish et aL]|2005| |Zapata et al. 2008 1 Walsh ^FaL] outflows associated or not with cm continuum (e.g., Argon et al. 



red-shifted emission, and the H2 emission. 



1998), the CO(2-l) HV blue-shifted emission, the SO blue- and 2003 1. In the case of OH masers tracing molecular outflows, the 



non detection of an ultra-compact Hn region can be explained 
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Fig. 6. Continuum emission at 1.3 mm observed with the SMA (Leurini et al., in prep.) from the maser zone in IRAS 17233 (solid 
black contours from 0.12 Jy beam 1 (3cr) in step of 0.24 Jy beam -1 ). Overlaid on the map are the positions of the OH (blue and 
Fish et aL 2005 1, of the CH3OH (black triangles, [Walsh et"aL 1998| l and of the water (blue and red open triangles, 
'2008j > rriasers. The green c ontour levels show the continuum emission at 1.3 cm (from 1.4 mJy beam 1 (3cr) in step 
of 1 .4 mJy beam" 1 , Zapata et al. 2008 1, while the grey scale map is the H2 emission. The dashed blue contour levels represent 



red squares, 



Zapata et al. 



the CO(2-l) blue-shifted emission at high velocities (from 5cr in step of 20cr, see Table [2]), the solid blue and red contours the 
SO(56 - 45) blue- and red-shifted emission (blue: from 80 to 104 Jy beam" 1 in step of 10 Jy beam" 1 ; red: 10 and 12 Jy beam" 1 ). 
The black square marks the Rl position, detected in the EHV red-shifted CO(2-l) emission (thik solid red contour, as in Fig. [3]>. 
The cross marks the position of the peak position of the 1.3 mm continuum emission. The dashed and dotted lines and the arrow are 
as in Fig. [3] OF1, OF2, and OF3 label the three identified outflows. 



with a very early evolutionary stage of the source, or with an 
intermediate mass object which will never produce a detectable 
ultra-compact Hn region. In our case, VLA 2d could be the driv- 
ing source of this outflow. Finally, we have a pattern formed by 
red- and blue-shifted FbO masers (empty triangles), which are 
well known tracers of molecular jets/outflows, spread around the 
1.3 mm continuum peak. 

Given the current resolution of our CO observations and the 
multiplicity of the high-velocity clumps, we cannot derive the 



exact number of molecular outflows in the region, in particular 
those extended along the N-S direction. However, using the high 
angular resolution H2 and maser emission data, we can try to 
draw a scenario where at least three molecular outflows originate 
from the maser zone in IRAS 17233: 

1. the first one (IRAS 17233-OF1, indicated by a solid arrow 
in Figs. [3] and [6]) is well traced by a H2 jet which seems to 
arise from the 1.3 mm continuum emission peak: the driving 
source could be either VLA 2a or VLA 2b. In addition, the 
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Fig. 7. Maps of the integrated blue- and red-shifted emission in the 13 CO(2-l), C ls O(2-l) and SO(56 - 4s) transitions. The solid 
contours show the blue-shifted emission in the velocity ranges v=[— 55, -20] km s _I for the I3 CO(2-l) line, v=[-20, -11] km s _1 
for C 18 0(2-l) and v=[-40, -11] km s for SO(56 - 4s), and the red-shifted emission in the velocity range v=[15, 20] km s _1 for 
the 13 CO(2-l) line, v=[5, 10] km s" 1 for C I8 0(2-1) and v=[15, 25] km s _1 for SO(5 6 - 4 5 ). The contours start from a 5cr level and 
are in step of 5<x. The grey scale map is the continuum emission at 1.3 mm (from 0.12 Jy beam -1 (3cr) in step of 0.24 Jy beam -1 ). 
The symbols are as in Fig. [3] The synthesised beam is shown in the bottom left corner. 



H2O spots are well in agreement with a molecular outflow 
directed along the NW(red)-SE(blue) direction. The present 
12 CO maps confirm this possibility, with high-velocity blue 
emission located SE of the 1.3 mm peak (Fig.|6|l, and the Rl 
red clump which is located ahead the H2 NW jet. This pic- 
ture is consistent with the scenario depicted by SO, a well- 
known molecule whose abundance is enhanced in molecular 
outflows (e.g., Bachiller & Perez Gutierrez] 1997| Codella 



Jet al.||2005) : blue-shifted emission east of the 1.3 mm peak, 
and red-shifted emission in the NW (Fig. |6j. Blue-shifted 
emission is also detected in 13 CO and C 18 (Fig. [7ji in the 
SE of the 1.3 mm peak. Thus, the Rl clump seems to trace 
EHV gas associated with the mass loss traced by H2. The 
scenario for the SE blue lobe is less clear: we have high ve- 
locity blue emission but no EHV clump and, surprisingly, 
there is no H2 counterpart, when usually the blue H2 emis- 
sion is brighter than the red one due to extinction. A possible 
explanation would be that the blue-shifted jet is leaving the 
high density clump hosting the YSOs, moving in a gas where 
the density is so low that Fb is not excited. A more realistic 
solution is that the driving source of the H2 NW jet (VLA 2a 
or VLA 2b) had formed inside the 1.3 mmm core, but close 
to its far edge; in this case the red-shifted jet would travel in a 



less dense ambient medium, while the blue-shifted jet would 
travel deeply into the molecular cloud with a consequently 
higher extinction. 

2. a second outflow can be found along the direction traced by 
the B3-B1-R2-R4 clumps (IRAS 17233-OF2, dashed fine in 
Figs. [3] and [6}. Note that this direction is well aligned with 
both the OH bipolar outflow and definitely supports VLA 2d 
as driving source; 

3. the third one can be found along the path traced by the B4- 
B2-R2-R4 clumps (IRAS 17233-OF3). In this case, the cor- 
responding dashed line of Figs. [3] and [6] passes through the 
H2O masers and close to the VLA 2a/1.3mm source, which 
could then host the driving source. 



This picture has to be considered as a first incomplete step 
to prepare higher angular resolution observations, in CO or in 
a classical tracer of the hot jet component as SiO, which are 
needed for a better understanding of the number of outflows in 
the region, and for the identification of the powering sources. 
For instance, there is no clear solution accounting for the R3 
clump. A possible blue counterpart of R3 could be found in B3 
and the driving source should be consequently in the 1.3 mm 
object. Similarly, R2 and R4 are not clearly assigned to any given 
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outflows, but in our current picture they could either belong to 
IRAS 17233-OF2 or to IRAS 17233-OF3. 

In conclusion, although higher angular resolution data are 
needed for a better understanding of the number of outflows in 
the region, and for the identification of the powering sources, 
on the basis of the existing data we can safely conclude that 
at least three molecular outflows are found in the maser zone 
in IRAS 17233. Given their association with CH 3 OH and OH 
maser emission, and given the presence of at least four compact 
sources of cm continuum emission, we conclude that the molec- 
ular outflows originate from intermediate or high-mass YSOs. 



5. Derived outflow parameters 

Global properties of outflows can be derived, under given as- 
sumptions, from the high-velocity CO emission. We will use the 
13 CO(2-l), 12 CO(2-l) and 12 CO(6-5) data to infer the optical 
depth and the excitation temperature. However, there is a major 
caveat in our analysis: interferometric data are affected by miss- 
ing flux from extended structures. Moreover, optically thick lines 
have a more uniform distribution than optically thin transitions, 
and may be more sensitive to interferometric spatial filtering. 
Thus, in the following discussion we may (1) overestimate the 
optical depth of the I3 CO(2-l) line; (2) overestimate the ratio 
between the 12 CO(6-5) and the 12 CO(2-l) lines, and therefore 
(3) overestimate the excitation temperature of the gas. In both 
cases, the problem comes from the 12 CO(2-l) transition, whose 
observed flux may be a lower limit to its true value. Since no sin- 
gle dish 12 CO(2-l) observations of IRAS 17233 are available, 
we cannot further investigate this possibility. In section § 5.1 
we will derive the optical depth of the 13 CO(2-l) and 12 CO(2 



1) lines; in § 5.2 we will investigate the excitation conditions 
of the gas along the outflow. Finally, in § 5.3 we will derive the 



properties of the outflow making use of the opacities derived in 
§ |5.1| and for a range of temperatures compatible with the results 
of the statistical-equilibrium calculations performed in section 
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Fig. 8. Overlay of the opacity of the 13 CO wings (t^; grey scale, 
levels are from 0.07 (light grey) to 0.12 (dark grey) in step of 
0.01) as averaged in the -55,-20 km s' 1 (red) and +15, +25 km 
s (blue) ranges, and the I3 CO integrated blue- and red-shifted 
outflow (see Fig. [7] for the sake of clarity only selected contour 
levels are drawn). The triangle marks the position of the CH3OH 
methanol masers (see Fig. [6). 



5.1. Optical depth 

By comparing the intensity of the I3 CO and I2 CO(2-l) inten- 
sities in the wings, it is possible to estimate the optical depths 
of both lines, assuming a [ 12 CO/ 13 CO] relative abundance of 77 
Wils on & Rood| |199 4). We estimated the optical depth of the 
J CO(2-l) line, T13, under the assumptions that the 13 CO(2-l) 
and 12 CO(2-l) lines have the same excitation temperature, and 
that the 12 CO(2-l) transition is optically thick. Since (i) the S/N 
ratio in the 13 CO(2-l) wings is sufficiently high only in the rel- 
atively low velocity outflow range, and (ii) self-absorption af- 
fects the 12 CO(2-l) line at velocities close to the ambient one, 
we derived the optical depths in the -55,-20 km s _1 (red) and 
+ 15, +25 km s _1 (blue) ranges. Figure [8] shows the distribution 
of the optical depth of the 13 CO(2-l) line averaged over these 
velocity ranges: the emission in the I3 CO(2-l) line results to 
be optically thin, with t 13 < 0.1-0.3. The I2 CO(2-l) optical 
depth, T12, was derived by simply scaling the T13 value by the 
assumed [ 12 CO/ 13 CO] ratio, thus resulting in an optical depth 
for the 12 CO(2-l) line, t 12 , between 10 and 30. 

5.2. LVG analysis 

In order to constrain the physical conditions associated with the 
CO outflow(s) by means of statistical-equilibrium calculations, 



we compared the I2 CO J = 6-5 and 2-1 lines {E u = 116.2 and 
16.6 K, respectively). For a proper comparison, the 12 CO(2-l) 
image, originally obtained with a 5'.' 5 x 2'.'1 resolution, was 
reconstructed at the lower resolution of the 12 CO(6-5) map 
(HPBW = 974). We then measured the ratio between the 
brightness temperature of the J = 6-5 and 2-1 lines. Figure [9] 
shows the distribution of with velocity with respect to the 
systemic one (Vlsr) as observed towards the blue- (upper panel) 
and red-shifted (lower panel) outflow lobes. Interestingly, Fig. 
[9] shows that, for each blue- and red-shifted lobe the higher the 
distance of the clump from the driving source(s) is, the higher 
is the R(,2 ratio, thus indicating higher excitation conditions. For 
instance, the R(,2 ratio calculated for the B4 clumps is very high: 
~ 10. However, the CO emission is associated with multiple 
outflows, and among the CO peaks only B2 and B4 clearly ap- 
pear to be related to the same flow. Therefore, the trend seen in 
the /?62 ratio ma y not be physically relevant for the other posi- 
tions. Higher angular resolution observations are needed to get a 
clearer picture of the region. In addition, Fig.|9]suggests also, at 
least for the red clumps, an increase of the excitation at the high- 
est velocities, probably closely associated with recently shocked 
gas. 

In order to constrain the excitation conditions, we used a 
Large Velocity Gradient (LVG) code with collisional rates from 
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Fig. 9. Distribution with velocity with respect to the systemic 
one (Vlsr) of the ratio between the brightness temperatures 
of the J - 6-5 and 2-1 lines as observed towards the blue- (upper 
panel) and red-shifted (lower panel) outflow lobes. 
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Fig. 10. Results of statistical equilibrium calculations for 12 CO. 
The ratio R(,2 (grey scale) between the brightness temperature of 
the J = 6-5 and 2-1 lines, as a function of the T^ n x «h 2 product, 
which can be used as a measure of the gas pressure, and 7\; n . We 
have assumed the LVG optical depth parameter Nco Av~' = 10 16 
cm -2 (km s -1 ) -1 . Solid lines mark R^ as measured towards the 
red- (upper panel) and blue-shifted (lower panel) outflow lobes. 
Dot-dashed lines are for the R^ uncertainties. Note that the R^ 
value towards B4 is out of the sampled plane (~ 13): white ar- 
rows are drawn to clarify the portion of the plane allowed for 
B4. 



Schoier et al. ( 2005 ). These include extrapolation to energy lev- 
els up to J=40 and collisional temperatures up to 2000 K, based 



on the original datasets calculated by Flower (2001 ) and Wernli 



et al. ( 2006| l for temperatures in the range from 5 to 400 K and 
energy levels up to J=29 and J=20 for collisions with para-H2 
and ortho-H2, respectively. 

We assumed a specific column density Nrg Av _1 ranging 



from 10 to 10 cm (km s ) . Figure 10 shows our typi 



conditions of the outflow clumps. Regarding the red lobe, the 
Rl and R2 clumps have 7\j n > 50 K and a product 7\j n x n Hl , 
which can be considered a measure of the gas pressure, > 10 6 K 
cm -3 . A similar scenario can be found in the blue lobes Bl and 
B2. On the other hand, the R^ value towards the extreme high 
velocity position B4 is out of the plane sampled by our models 
(,^62 ~ 10). The red-shifted extreme high velocity position (R3, 
R(,2 ~ 7) is at the limit of the plane. Such high values of R(,2 sug- 
gest very high temperatures (> 150 K for R3 and > 300 K for 
B4) for the gas at those positions; however, these values must be 
taken as upper limits since the flux of the CO(2-l) line could be 
underestimated due to filtering of large structures, as stated in 
Sect. 5. 



cal model results, illustrated graphically with the ratio R^ (grey 
scale) as a function of the thermal pressure (7\i n x «h,) and 
the kinematic temperature (7\i„). Figure [10] shows the solu- 
tions found for A^co Av~' = 10 16 cnT 2 (km s -1 ) -1 : however, 
this should be considered as an upper limit since also solu- 

tions down to 10" cm 2 (km s 1 ) 1 are acceptable. Solid lines &a 0utffow e ef/cs and kinematics 
mark the measured values of R^ towards the red- (upper panel) 
and blue-shifted (lower panel) outflow lobes. For the sake of 
clarity we plotted the values averaged over the whole velocity 
ranges. In addition, we further averaged the measurements of the 
clumps associated with the lower excitation conditions, showing 
similar values (red: fl 62 (Rl,R2) = (R 62 (Rl)+R(, 2 (R2))/2; blue: 
fl 62 (Bl,B2) = (fl 6 2(Bl)+7? 62 (B2))/2, where R 62 at each posi- 
tion is the average over the values for different velocities). Dot- 
dashed lines are for the R^ uncertainties, defined as the mini- 
mum and maximum values that /?62(R1,R2) reaches. 



Figure 10 clearly shows that we cannot obtain severe con- 
straints on the physical conditions of the gas with only two tran- 
sitions. However, it allows to quantify the different excitation 



The outflow mass, moment and energy (M out , p out , and E out , re- 
spectively) were derived from line emission using the classical 
formulae (e.g., |Lada| [T985] ) as a sum over all the velocity chan- 
nels within the blue and red wings and within the areas defined 
by the 3cr contour levels of the respective outflow lobes. We 
sum all the individual pixel contributions within the lobe area. 
In other words, the sum of the contributions from each pixel 
and velocity channel builds up the total outflow mass, moment 
and energy. We used the 12 CO(2-l), I2 CO(6-5), and I3 CO(2-l) 
lines separately in order to have different estimates. The abun- 
dance ratio r 12 CQl/rH ? l and ["COl/ TH?! was taken as 10~ 4 and 
1 .3x 10~ 6 (e.g., Scoville et al. 1986 1, respectively. We conserva- 
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Table 3. Whole emission (summing the contribution due to the 
three identified outflows IRAS 17233-OF1, -OF2 and -OF3) 
and the IRAS 17233-OF1 outflow energetics, derived assuming 
distance of 1000 pc and an excitation temperature in the 50-200 
K range (see text) 



Tracer 


T/n 1 ^out* Pout ' 


F * 

^OUt 




(M e ) (Mokms- 1 ) 


(10 45 erg) 


Whole emission (OF1+OF2+OF3) 


12 CO(2-l) 


19.0 2.4-3.7 99-140 


60-90 


13 CO(2-l) 


0.3 1.0-1.5 19-30 


4-6 


12 CO(6-5) 


0.1 0.7-1.1 20-29 


6-10 


OF1 red-shifted outflow 


12 CO(2-l) 


23.0 0.6-0.9 18-28 


6-9 



Whole emission: mean r between blue and red emission optical 
depths (Sect. 5.1) 

Values derived for T cx 50 and 200 K (Sect. 5.2) 



tively assumed an excitation temperature T ex in the range 50-200 
K for all velocities, following the LVG analysis. We corrected 
the M out , p 0M , and E out for the optical depths derived in Sect. 
5.1. The average distance of 1000 pc has been used. The 12 CO 
molecule is the usual outflow tracer due to its high abundance, 
which allows to detect the weakest wings at the highest velocity. 
Nevertheless, the risk of 12 CO lies in its high opacity, which can 
bias the observations towards the colder and more extended por- 
tion of the outflow. We tried to verify this possibility by using 
also the 12 CO(6-5) and 13 CO(2-l) lines, which present weaker 
emission, thus preventing the detection of the highest velocities, 
but which could possibly better trace the more collimated out- 
flow component. 

Table [3] lists both the values derived summing the contribu- 
tion due to all the outflows present in the region (hereafter whole 
emission) as well as the values of the OF1 red-shifted outflow. 
In the latter case, we refer to the 12 CO(2-l) line only, since it 
is the only transition where the OF1 outflow is clearly detected. 
Note that, for the whole emission, we do not indicate Tbi ue and 
T re d, but, for the sake of clarity, the mean value, r mean = (Tbiue + 
T re d)/2, which can differ slightly. However, the average values of 
the optical depth in the red- and blue-shifted lobes, Tbi ue and r re d, 
were used in the calculations. For the 12 CO(2-l) line, we adopt 
Tbiue = 15 and r re d = 23; for the 13 CO(2-l) transition, Tbiue = 0.2 
and T re d = 0.3; for the 12 CO(6-5) line, we used a typical value 
inferred by the LVG analysis, Tbi ue = T re d = 0.1 

From Table [3] we confirm what expected, i.e. that the 
12 CO(2-l) line allows us to trace the most energetic portion of 
the outflow(s) thanks to the detection of the highest velocities. 
The I2 CO(6-5) and 13 CO(2-l) lines report E out lower by more 
than one order of magnitude. Note also that 13 CO is tracing a 
mass ~ 2-3 times lower than that derived by the 12 CO(2-l) tran- 
sition. This is consistent with the fact that only limited velocity 
ranges of 13 CO(2-l) can be used for the integrated intensity in 
the red-shifted emission due to contamination from several HC 
tracers (see Fig. [2 1, and with the fact that the 13 CO emission 
is more compact than that of the I2 CO(2-l) line. In addition, 
we cannot exclude that the mass estimate is affected by other 
uncertainties due to the presence of temperature gradients, with 
the 13 CO tracing warmer inner regions, as well as to the opti- 
cal depth estimate. Therefore, we derived the outflow kinematics 
only by using the I2 CO(2-l) emission. 

In order to derive the real outflow parameters, the angle 
8 of inclination to the plane of the sky has to be assumed. 



Unfortunately, because of the complexity of the region and the 
outflow multiplicity, we cannot derive a good estimate of 6. 
However, taking into account the quite elongated structure of 
the outflow(s) and the very limited overlap between red and blue 
lobes, it is reasonable to assume an intermediate inclination, in 
the 30° -60° range. A kinematic age for each flow lobe, f^n, was 
determined from the projected distance between the position of 
the outflow lobe and that of the mm-peak, and the correspond- 
ing radial outflow velocity component. From the kinematic out- 
flow age, the mass entrainment rate of the molecular outflow, 
M oul , the momentum rate, p, ot , and the mechanical luminosity, 
L mec , were derived for the blue and red lobes separately, and 
then added together to obtain the total values presented in Table 
[4] There, we list fkj n , M out , p tot , L mec without correction for the 
inclination angle, but also for equals 30° and 60°. As done 
for optical depths, we give a mean value for the kinematic age 
between the values derived separately for the blue and red emis- 
sion. After geometry correction, for the whole emission we have 
a kinematic age of 10 2 -10 3 yr, and consequently M out ^ 2-9 
10~ 3 M Q yr 1 , p tQt 0.1-0.7 M km s" 1 yr 1 , and L mec 500- 
8000 Lq. By comparing these estimates with the expected cor- 
relation between the dynamical parameters and the source bolo- 
metric luminosity (e.g., Beuther et al. 2002b| l, and taking into 
account the outflow multiplicity, we infer for the driving sources 
a Lboi ^ 10 4 L , i.e. earlier than a B0. 5 type. On the other hand, 
if we take into account the OF1 outflow only, the bolometric lu- 
minosity of the driving source is of the order of 10 4 L G . The 
luminosity estimated for the driving sources are well in agree- 
ment with the bolometric luminosity of IRAS 17233, which is 
2.5 xlO 4 L for a distance of 1 kpc. 

The derived values are in general good agreement with those 
(M out = 0.8 M , p out = 25.3 Mq km s" 1 , £ out = 8 .6 x 10 45 erg, 
fidn = 6200 yr) reported by |Leurini et al. ((2008) based on the 
analysis of single dish data in CO(3-2), and with the conclu- 
sion that the kinematical outflow parameters are typical of mas- 
sive YSOs. The differences between the two analyses can be ex- 
plained with the different areas and different velocities adopted 
for the analysis. Moreover, our previous study on CO(3-2) did 
not include any correction for the opacity of the line. 

5.4. Age vs. outflow collimation 

The collimation and the age of the outflows can be compared 
with the scenario described by Beuther & Shepherd (2005) 
where jet-like outflows occur only in early evolutionary phases 
where no HCHn region has formed and with stellar luminosi- 
ties corresponding to late B-type main sequence values, whereas 
in later evolutionary phases, a HCHn forms and the wind from 
the central massive star produces an additional less collimated 
outflow component. Since to form early O-type massive stars 
via accretion the protostellar objects must accrete even after the 
central object has reached the main sequence, the evolutionary 
scenario depicted above also corresponds to a change of lumi- 
nosity of the (proto)star during its evolution, from a late B-type 
to an early O-type star before reaching its final mass and stellar 
luminosity. This scenario would explain the observational result 
that no collimated jet-like outflow has ever been detected from 
very young early O-type (proto)stars. 

In the IRAS 17233 case, the N-S outflows show a remarkable 
degree of collimation ( f c ~ 4) as well as a relatively small open- 
ing angle, 5 ~ 11°. In addition, these estimates have to be con- 
sidered as lower (collimation) and upper (opening angle) limits, 
given the unknown outflow multiplicity. The spectral type (B0. 5) 
derived by the outflow parameters and the age of flow (~ 10 2 - 
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Table 4. Whole emission (due to the three identified outflows 
IRAS 17233-OF1, -OF2 and -OF3) and IRAS 17233-OF1 out- 
flow timescales and kinematics, derived from 12 CO(2-l) emis- 
sion, assuming distance of 1000 pc (see text). Uncorrected val- 
ues as well as values corrected for assumed inclinations to the 
plane of the sky (9) are reported 



e 




M J p t ot* 




(deg) 


(yr) 


(M Q yr" 1 ) (M Q km s" 1 yr" 1 ) 


(£o) 


Whole emission (OF1+OF2+OF3) 




480-940 


3-5 10- J 0.1-0.2 


760-1200 


30 


278-545 


5-9 10~ 3 0.4-0.7 


5244-8280 


60 


830-1626 


2-3 10~ 3 0.09-0.13 


585-924 


OF1 red-shifted outflow 




360-600 


1-2 10" 3 0.03-0.05 


86-132 


30 


209-348 


2-3 10~ 3 0.1-0.2 


593-911 


60 


623-1038 


~ 1 10~ 3 0.02-0.03 


66-102 



1 Whole emission: mean /kin between blue and red emission (Sect. 
5.1) 

* Values derived for T cx 50 and 200 K (Sect. 5.2) 



10 3 yr) agree with the scenario of Be uther & Shepherd| (2005). 
However, the OH maser spots trace a bipolar outflow with a def- 
initely larger opening angle (see Fig. [6]i than that observed in 
CO and H2, thus suggesting we are observing a stage where jet- 
like and less collimated winds co-exist. The detection of HCH11, 
typical of more evolved objects, also supports our interpretation 
that IRAS 17233 is in a transitional phase, when the star also 
produces an additional less collimated outflow component. 

6. Conclusions 

We presented interferometric (SMA) and single-dish (APEX) 
observations of one of the nearest massive star forming region 
IRAS 17233 in CO(2-l) and (6-5) with angular resolution of 
5'.'4 x 179 and 974, respectively. The main results can be sum- 
marised as follows: 

1. The data reveal a clumpy extended (~ 50") structure with 
well separated blue- and red-shifted emission, and an overall 
structure roughly aligned along the N-S direction and cen- 
tred on the region where groups of H2O, CH3OH, and OH 
maser spots are detected. This region is also associated with 
the peak of the 1.3 mm emission we mapped, as well as, 
on subarcsec-scale, with 4 HCHn regions (called VLA 2a, 
b, c and d) observed at cm-wavelengths and indicating mas- 
sive star formation. The outflow is associated with extremely 
high velocities, up to — 200 and +120 km s _I with respect to 
the ambient LSR velocity. Up to eight outflow clumps have 
been observed in CO(2-l). The whole outflow structure is 
very well collimated ( f c ~ 4) and with an opening angle of 
-11°. 

2. The large number of outflow clumps clearly indicates out- 
flow multiplicity, as expected for a cluster of high-mass 
YSOs. Using high-angular resolution maps of H2, cm- 
continuum, and maser emission, we tried to identify different 
outflow components. At least three outflows can be individu- 
ated. One outflow (IRAS 17233-OF1) is compact (~ 5"-10") 
and located along the NW-SE direction: it shows extremely 
high-velocity red-shifted emission and it is associated with 
a jet traced by H2 emission as well as with red and blue 
H2O maser spots and SO emission. Its driving source has 



to be found in the region associated with the 1.3 mm ob- 
ject, the CH3OH masers, and two massive YSOs (VLA 2a 
and VLA 2b). On the other hand, the extended (> 50"; > 
0.2 pc) emission along the N-S direction can be disentangled 
into at least two more bipolar outflows: one (IRAS 17233- 
OF2) associated with a counterpart on smaller scales (~ 5"; 
~5000 AU) traced by red- and blue-shifted OH masers and 
driven by a massive YSO traced by cm-emission (VLA 2c or 
VLA 2d), and another one (IRAS 17233-OF3) whose driving 
source has to be found again among VLA 2a and the 1.3 mm 
peak. 

H2 emission is found mostly in correspondence to the blue- 
shifted clumps, probably due to extinction hiding the red- 
shifted H2 lobes. In two cases, we found extremely high ve- 
locity CO emission associated with H2 emission, suggesting 
we are observing material accelerated along the jet traced by 
H 2 . 

The molecular outflows were observed also in isotopologues 
of 12 CO: i.e. 13 CO(2-l) and C I8 0(2-1), as well as in SO(6 5 - 
54). Although covering a smaller velocity range with respect 
to the 12 CO(2-l) line (due to lower line brightness), also in 
this case the outflow structure is well traced. By comparing 
the intensity of the 13 CO and I2 CO(2-l) wings, we have es- 
timated the optical depths: the emission in the I3 CO(2-l) 
line results optically thin, with T13 < 0.1-0.3. Consequently, 
the 12 CO(2-l) optical depth (derived simply scaling by the 
assumed [ 12 CO/ 13 CO] ratio) is quite high (-10-30), indicat- 
ing thick emission even at outflow velocities + 30-60 km s" 1 
with respect to the ambient emission. 

To estimate the relative 12 CO excitation conditions of the 
outflow clumps, we compared the J = 6-5 and 2-1 line 
brightnesses. For the blue-shifted lobe the higher the dis- 
tance of the clump from the driving source(s), the higher is 
the /?62 ratio, suggesting higher excitation conditions. In ad- 
dition, at least for the red clumps, we find indication of an 
increase of the excitation at the highest velocities, probably 
closely associated with recently shocked gas. 
Although only with two transitions, LVG statistical- 
equilibrium calculations give a rough estimate of the physi- 
cal parameters of the molecular gas traced by CO. The red 
and blue clumps closer to the YSOs are associated with 7\i n 
> 50 K and with a gas pressure > 10 6 K cirT 3 . 
The estimate of the kinematical outflow parameters such as 
momentum, kinetic energy, and the mechanical luminosity 
show values that, even taking into account the outflow mul- 
tiplicity, are typical of massive YSOs (Lb i > 10 4 L Q ), and in 
agreement with the measured whole bolometric luminosity 
of the source (Lboi = 2.5 10 4 L Q ). 

The kinematic ages of the flows are in the range 10 2 - 10 3 yr, 
and therefore point to young objects that still did not reach 
the main sequence. We compared our results with the sce- 
nario described by Beuther & Shepherd (2005) where jet- 
like outflows are powered only by young protostars with lu- 
minosities that correspond to late B-type main sequence val- 
ues, and still without an HCHn region. In the case of the N-S 
outflows in IRAS 17233, we observe a high degree of colli- 
mation as well as the indication of definitely less collimated 
wind traced by OH maser spots. We thus suggest that we are 
observing a stage where jet-like and less collimated wind co- 
exist, in agreement with the detection of very small (~ 10~ 3 
pc) ionised regions and by the spectral type (B0. 5) derived 
from the outflow parameters. 

Our analysis shows the importance of high sensitivity obser- 
vations to detect extremely high velocity emission to access 



14 



S. Leurini et al.: Extremely high velocity gas from the massive YSOs in IRAS 17233-3606 



the most energetic part of molecular outflows and detect the 
counterpart of the H2 emission. Moreover, high-spatial reso- 
lution observations of higher excitation CO lines as well as 
SiO, sensitive to the hot jet component and not contaminated 
by ambient gas emission as low rotational CO transitions, 
would be instructive to better assess the outflow multiplicity 
in the source and for a better characterisation of the excita- 
tion conditions. 
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